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O 

We show that it is possible to obtain numerical solutions to quantum mechanical problems in- 
volving a fractional Laplacian, using a collocation approach based on Little Sine Functions (LSF), 
, which discretizes the Schrodinger equation on a uniform grid. The different boundary conditions 

1) ■ are naturally implemented using sets of functions with the appropriate behavior. Good convergence 

properties are observed. A comparison with results based on a WKB analysis is performed. 
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PACS numbers: 03.65.Gc,02.70.Jn,11.15.Tk 



I. INTRODUCTION 



C$ ' There has recently been great interest in what is called fractional quantum mechanics. Laskin [l|, 0, H[ derived a 
fractional Schrodinger equation from a fractional version of the path integral. The nature of the fractional quantum 
mechanics is determined by the Levy index < a < 2 and the requirement for the first moment's existence gives the 
restriction 1 < a < 2 Q. The author first solved the fractional Schrodinger equation for the infinite potential well, 
the Bohr atom and introduced a fractional oscillator Q. Later he derived the fractional Schrodinger equation for 
^ | three-dimensional motion and solved it for the Bohr atom and a one-dimensional oscillator by means of semiclassical 

■ approaches 0]. 

Guo and Xu discussed the solutions of the fractional Schrodinger equation for the free particle, the infinite well, 
, and a simple model for barrier penetration, among other physical applications Q and Dong and Xu [|| solved some 
• ' other examples in the momentum representation. 

Zoia et al Q addressed the problem of the boundary conditions in fractional Laplacian equations and proposed a 

■ method for the accurate calculation of eigenvalues and eigenfunctions that overcomes the difficulty that had arisen 
\ in earlier approaches when a — ► 2. In particular, they considered absorbing and free boundary conditions and took 

• • ■ advantage of the fact that the Laplacian equations exhibit exact solutions when a is an even integer in order to test 
their method. 

The purpose of this paper is to show that a simple collocation method is suitable for the treatment of the fractional 
' Schrodinger equation. This approach has already proved successful in standard quantum mechanics as well as for 
several other physical problems [1, 0, [13, HH, 03 ■ We think that it is most convenient to have a method that applies 
to a wide variety of problems and for this reason, in this paper, we propose the ubiquitous collocation method based 
on little sine functions (LSF) that easily accommodate to a variety of boundary conditions [|[. 

In Section |TT] we outline the main features of fractional quantum mechanics and discuss the fractional Laplacian 
operator in terms of well known operator methods. In Section IIIII we describe the collocation method, develop the 
Fourier decomposition of the sampling functions to be used, as well as calculate explicitly the effect of the fractional 
differentiation operators on such sampling functions. In section ITVl we apply our collocation method to two fractional 
differentiation problems. Finally, in section [V] we present our conclusions. 
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II. FRACTIONAL QUANTUM MECHANICS 

By means of the fractional path integral in which the Levy motion substitutes the Brownian one, Laskin 

0,1 

derived the fractional Schrodingcr equation 

ih^ = H a i), H a = -D a (HV) a + V(x), (1) 

where D a is a generalized fractional diffusion coefficient. The author proved that the fractional Schrodingcr operator 
H a is Hermitian or self-adjoint 0, Si]. The meaning of the fractional derivative is clear from its effect upon a plane 
wave 0,0,3 

\7 a e ikx = \k\ a e lkx . (2) 
Laskin Q also considered a 3D generalization of the fractional derivative: 

-V Q ^(-A) Q / 2 , (3) 

where A is the Laplacian operator. 

On the other hand, Zoia et al [1] considered the alternative definition 

_^_ p iqx = _|g|« e <9* ] 



d\x 



= -(-A)"/ 2 , (4) 



d\x 

that is consistent with the one above. Clearly, the effect of the fractional derivative on an arbitrary function follows 
from the form of the standard Fourier transform 0, 0, H, 0, HI, 0] ■ 

In this paper we resort to the standard definition of the function of an operator in order to define the fractional 
derivative. One advantage of this definition, which is consistent with the one just outlined above, is that it enables 
us to obtain several results without further proof. If A is a Hermitian operator with a complete set of eigenvectors 
{\J >}, 

A\j>=aj\j>, (5) 
then we define the function f(A) by means of the spectral decomposition [l3[ 

/(i)=^/(a i )|i><J/1, (6) 

3 

that makes sense if f(x) is well defined for all x = a,j. The extension to a continuous spectrum is straightforward. 
Notice that if f(x) is real then f(A) is Hermitian. Besides, if [A, B] = then, obviously, [f(A), B] = 0. 

For example, on the whole line we have the standard definition in terms of the Fourier transform 0, 0, 0, 0, M, 

1 />oc />oo 

(-Ar^(x)^-I I \k\ a e lk ^dk^{y)dy, (7) 



— oo J — oo 



and, in terms of the Fourier series for periodic boundary conditions, 

(~A) a /^(x) = — J2 \<z m[x - v H{y) dy. (8) 



n— — oo 



In Section lllll wc consider other functions and other boundary conditions. 

In addition to being the basis for the numerical method proposed in this paper, the well known general operator 
results outlined above enable us to derive several conclusions without further proof. For example, if we choose the 
domain of the operator — A so that it is Hermitian, then its eigenvalues are real and positive. Its obvious consequence 
is that (-A)"/ 2 is also Hermitian 0,0]. 

It is well known that the Laplacian commutes with the operator that produces the inversion transformation r — * — r; 
therefore, (— A) Q / 2 also commutes with the inversion operator Q. 
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Every textbook on quantum mechanics shows that the nontrivial solutions to the eigenvalue equation (— Aip)(x) 
-ip"(x) = Xip(x) with boundary conditions ip(±a) = are given by 



A = A 



1 



: sin 



nir(x + a) 



Aa 2 



2n 

; 71= 1,2 



(9) 



Therefore, it follows from the operator equations outlined above that 

f a ip n = E%ip n , 

f a = D a (-h 2 A) a/2 , E" = D a 



~2a 



(10) 



In other words, it is not necessary to solve the time-independent fractional Schrodinger equation explicitly 0, H| 
because we already know that the eigenfunctions are exactly those of the particle in an infinite well. 
As a final example, consider the time-dependent Schrodinger equation 



(11) 



with the boundary conditions ip(±a, t) = 0. Straightforward application of the method of separation of variables leads 
to 



n=l J - a 



(x)*ip(x, 0) dx, 



(12) 



that generalizes the result derived earlier for the same fractional model 0, H| . 

Notice that present definition of the fractional derivative avoids the problem of nonlocality of the Riesz derivative 
pointed out by Jeng et al [Til ]. 



III. THE METHOD 



As stated above, our purpose is the application of a collocation method to fractional differentiation problems on 
bounded intervals, with a variety of boundary conditions fll| . 

This collocation method starts by sampling a function / on [— L,L] in a finite uniform set of points Xk, and then 
interpolating by 

N 

J2 f(xk)s k (N,L,x), 

k=-N 

where the sampling functions Sk{N,L,x), called little sine functions (LSF). They arc defined as 

/at- r s S N (x,x k ) 

s k {N,L,x) = — -, 

d N (xk,Xk) 

where 5jy(x,y) — X^=o 4>n(x)4>n{y) is the partial sum kernel for a complete orthonormal set of functions {<fi n } in 
L 2 ([—L,L]), suitable chosen according to given boundary conditions [ill ]. 
Thus, for a given operator T, we define Tf on [-L, L] by 

N 

Tf(x)= Y, f(xk)Ts k (N,L,x), 

k=-N 

which then can be calculated by interpolating, in turn, the function Tsk{N, L,x). As we are interested in the 
application of T when given as a function of differentiation, it is convenient to decompose each sampling function 
Sk(N, L, x) in exponential functions. 

We therefore develop the Fourier decomposition for the sampling functions Sk{N, L, x), taken from the LSF sets in 
[Tl| . in order to apply fractional differentiation. 
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A. Periodic boundary conditions 

When given periodic boundary conditions, f(—L) = f(L), we use as an orthonormal set the functions 

2Lk 

with sampling points x k = 7^ — -, \k\ < N. We will denote the corresponding set of sampling functions s k (N, L, x) 
as LSFi. 

In order to obtain the Fourier decomposition for s k (N, L,x), we start by calculating 

N N 

One can verify that, for any y G [— £, L], Sjy(y, y) = — —J- — and hence 

N 



r(a; ' y) = 2ivTT £ 1 



2£fc 

Therefore, since x k = ^ ^ ^ , |fe| < AT, 



1 * 

s fc (iV, L, g) = 6 N (x, x k ) = 2N + 1 E 



e 2«+i e L 



Note that we can write 



1 W 2 *, 2 1 2N 1 + (-1)" h 

s k (N,L,x) = e ~ l ^TTe l ^ L:E = '—e~ l ^rfie l ^ : 

y ' ' ' 2N + 1 ^ 2N+1 ^ 2 

n=-N n=-2N 
2N 

AN + 2 ^ K v ; 

n=-2JV 

We prefer the use of the last expression, as it resembles the formulas below. 



B. Dirichlet boundary conditions 

In the case of given Dirichlet boundary conditions, f(—L) = f(L) = 0, we use the complete orthonormal system 



M*) = j= cos ( (2n ^ lTO ) , M*) = ^ "in ( JL 3 JE ) . 
Lk 



with sampling points x k = -j^-, \k\ < N. The corresponding set of sampling functions will be denoted by LSF2. 
Now, we have 

N-l ^ 2N 

5 N (x,x k ) = E {ipn(x)tp„(x k ) + <t>„{x)(f) n {x k )) = - ^ (e^T^r _ (_l)V w) e *ff ^ 

n=0 n=-2N 
TV 

and note that S n [x k ,x k ) = — . Thus 

L 



s k (N,L,x) = ^- E (e^w _ (_l)Vw) e iff= 



n=-2Af 
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We observe that 

e -*w — (_i) n e lJ Sr = e 1 ^ (e'^e -4 ^ — e^e 1 ^) = 2i n ~ 1 sin ^— + 

so we obtain 

2iV 



7rn 



1 fc 

2 + 2N 



C. Antiper iodic boundary conditions 

When given antipcriodic boundary conditions, /(— L) = —f(L), we take as a complete orthonormal system 
, f \ 1 /2n+l x 1 . /2n + l \ 



Lk 

with sampling points x k = — , |fc| < TV. We will denote the corresponding set of sampling functions by LSF 3 . 
This time we have 



S N (x,y) = J2 (Mx)Mv) + Mx)Mv)) =^7 ^ e l ^< x 

1 2N 



ti=0 n=-N 
2N 



--2N 



N 

since we are adding only on odd numbers. As 5n(u, y) = — for every y € [— L, L], we have 

L 



2N 

s k (N,L,x) = — J2 I 1 "("I)' 



47V 

n=-2Af 



D. Neumann boundary conditions 



In the case of given Neumann boundary conditions, f'(—L) = f'(L) = 0, we consider the complete orthonormal 
system 

= V2T' = VI C ° S \L /' "- = 1,2,..., 



= cos ( 2? 2^ 1 7TX 



with sampling points x k = -^rz — -, |fc| < N. We denote the corresponding sampling function set by LSF4. 
This time we have 



N 



S N (x,X k ) = 2J (^n(x)lpn(xk) + <Pn (x)<p„ {x k )) 
1 2N 

— ^2 (e^ 1 ^^ + (—l) n e l ^+^)e l ^ x , 



4L 

n=-2N 



and 5N(xk, %k) = — ^r- — ■ T 



2L 



j 2N 

s k (N,L,x) = 2 , 2N+l s J2 (e"''^ 1 + (-l)"e l ^TT)e^ ; 



n=-2N 



Using the identity 



e ^n'+i _(_ (_i)™ e I 2N' r +i — e »T ( e * T a r e 1 2n+i _|_ e l ™2" e'2»+i ) = 2i n cos 



1 



2 27V + 1 



we have 



s fe (7V,L,a;) 



1 



2iV 



27V- 



2 COS 



1 



2 27V +1 



=-2JV 

As a summary, we note that we have written all the sampling functions in the form 

2N 

s k (N,L,x)= J2 C n {k,N)e^ x , 



where the coefficients C n (k,N) are given by 



-2JV 



2(27V+ 1) 
{ 



C n {k,N) = { 



1- -1 " -i-fct 

; — e 2jv 

47V 

i n cos( 2 + 5iv+T) n7r 



LSF 2 
LSF 3 
LSF 4 . 



27V + 1 

Note that these coefficients do not depend on the length of the interval [-L, L\. 



E. Differential operators 

We now note that, after taking the derivative of each s k , we obtain 

, 2N 

^s k (N,L,x) = Yl C n{k,N)(i^)< 

n=-2N 

. d 

so, if p is the momentum operator p = —i — , we obtain 

ax 



ps k (N,L,x)= ^ C„(fc,7V)( 



2AT 

-mr\ im , 
— )e 2L 
2LJ 

-2N 



We then define, for a given function m, the operator m(p) on the functions s k by 

2JV 

m(p)8 k (N,L,x) = ^ Cn (k,N)m(^y^ x , 

n=-2N 

i.e., we define m(p) through the spectrum of p. 

We are particularly interested in the case of the fractional operator (—A)"' 2 , where A is the Laplacian A 
-p 2 . Thus, we have (-A)"/ 2 = |p| Q , and 

2N 

CJk. N\ f- 

2L 



2N | | 

{-A) a '*s k (N t L,x) = Yl C n (k,N)(^Ye^ 



In order to interpolate the resulting functions (— A) Q / 2 Sfc, we calculate its value in the sampling points Xj, i.e. 



2JV 



(-Ar^s k (N,L, Xj )= J2 C n (k,N)(^ 



\n\n\ a 



--2N 



We do this explicitly for each of the sampling function sets from above. 



(LSFi) We first calculate it for the set LSFi. In this case Xj 



2Lj 
2N + 1 



and 



and thus 



C, 



\n\ir\ a 1 + (-!)" i^k-j) 



2(2N+1) 



g 2iV + l 



so 



2N 



-A)<*/ 2 s k (N,L,x j ) = J2 



-2N 



|n|7T\« 1 + (-l) n _M£=i> 



2L I 2{2N + 1) 



g 2N + 1 



N 



2N 



n=l 



mt\ a 2n(k — j)n 



L 



2N+1 



(LSF 2 ) For this set, Xj = ^ and 



C n (k,N) 



2N 



1 k 

2 + 2N 



nix, 



so 







k \ i 


) 2^ Sm 













Thus we obtain 



2jV 



(-A) a / 2 8 k (N,L, Xj ) = J2 



nTT\ a i" 



=-2JV 
2N 

E 



2L J 2Ni 



sin 







[G- 


"2iv"H 







2N ^ \ 2L 

n—1 



/ n(k-j)n n(k + j)n 

COS rrr ( — 1) cos 



27V 



2iV 



(LSF3) For this set, Xj = and 



C n (A,JV) 



-1)" 



4iV 



e 2N 5 



so 



|n|7T\«l - (-1)™ ^(fc-j)* 

e 2JV 



2L 



47V 



and therefore 



2iV 



n=-2N 
N 



\n\n\ a l- (-1)™ inc».-j), 
e aw 



4JV 



(2n- 1 ) 7r A Q cos ( 2n ~ l)(fc-i)7r 



21 



2iV 
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(LSF4) Finally, for the set LSF4, the sampling points are given by Xj = 

i n (\ 

CJk, TV) = cos - 

" V ' 27V+1 \2 



2Lj 
27V + 1 

k 



2N+1 



and 



mr, 



so we have 



Thus 



C, 



2L ) 2N + 1 



cos 



1 



2 27V + 1 



2N 1 1 



r/i 



=-2JV 



27V + 1 



■ cos 



2 27V" +1 



g 27V+1 _ 



g 2JV + 1 



1 fnn\ (X ( n(k — i)ir 



— 1 ) cos 



n(k + j)7T 
27V + 1 



F. Collocation 

As described above, we approximate a function / on [-L, L] through interpolation from the sampling points x k by 
means of 



/(*)« E f{*k)sk{N,L,: 



k=-N 



Thus the action of m(p) is defined by 



m 



N 

(p)f{x)= f(x k )m(p)s k (N,L,x) 

k=-N 

N N 

= E f( Xk "> E m (P) s k(N,L,x :j )s : j(N,L,x) 

k=-N j=-iV 

TV , TV 27V s 

= E ( E E C n (A:,7V)m(^) e ^) Sj (7V,L,x). 

3=-iV ^k=-N n=-2N ' 

Note that we can view this as the action of the matrix 

M N p = [M N p] k , , —TV <k,j < TV, 
on the vectors f(xk), —TV < k < TV, where the matrix entries are given by 

2N 



[M N p] kj = ]T C n (k,N)m(^), 



= -2N 



Note that we have explicitly calculated these coefficients for (—A)"/ 2 = \p\ a in the previous section. 



IV. APPLICATIONS 



In this section we apply the present collocation method to the fractional versions of the anharmonic oscillator and 
the Mathieu equation. 
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A. Fractional oscillators 



Our first example is the fractional oscillator 

H = D a (-h 2 A) a / 2 + q 2 \r\P, (13) 

studied by Laskin [l[ by means of the semiclassical WKB approach. In the one-dimensional case Laskin obtained the 
following approximate analytical expression for the energies: 

En= { 2B(l/p,l/a+l) ) ( n+ 2) ' (14) 

where B(l/f3, l/a+ 1) is the beta function. 

This problem is suitable for illustrating the application of the collocation method described above. For concreteness 
we choose a dimensionless model with D a = q = h = 1 and resort to the set LSF2 because the Dirichlet boundary 
conditions are suitable for this problem. We obtain reasonably accurate results with N = 50, that corresponds to 99 
sampling points. 

The LSF are defined on an interval |x| < L. The unphysical parameter L may lead to inaccurate results if it 
is not chosen properly. If it is too small the wave function will decay too rapidly. If, on the other hand, L is too 
large then we would need an unnecessarily large number N of sampling points in order to have sufficiently accurate 
results. In order to get a reasonable balance between those parameters of the LSF method we resort to the strategy 
followed in earlier applications of the collocation approach based on sine functions [l5| and then extended to LSF (8[ . 
Since the trace of the Hamiltonian matrix is invariant under unitary transformations, and the actual eigenvalues are 
independent of L, then it is reasonable to choose the value of L close to a stationary point. This principle of minimal 
sensitivity (PMS) [l6[ gives the optimal value of L for a given value of N. In the present case the stationary point is 
a minimum that we will call Lpms from now on. 

In Table Q] we report the three lowest eigenvalues of the fractional harmonic oscillator ((3 = 2) with a = 3/2 
for grids of varying size. We appreciate that the value of Lpms shown in the second column grows with ./V in 
agreement with the argument above. The last row shows the energies obtained by application of exactly the same 
approach to the Schrodinger equation in the momentum representation, where no fractional derivatives are present, 
for a much finer grid. Notice that the rate of convergence of the present collocation method depends on a. In order 
to appreciate this point more clearly we may compare present results for N = 100, which exhibit only 4 exact digits, 
with identical calculation for the standard harmonic oscillator (a = 2), where just N — 10 enables us to obtain 
E PMS w 0.9999999999991 that is about 10" 13 off the exact value E% xact = 1. 



TABLE I: Three lowest energy eigenvalues of the fractional harmonic oscillator (f3 — 2) with a = 3/2. The last row shows the 
results for the Schrodinger equation in momentum space. 



N 


LpMS 


E 


E 1 


E 2 


10 


4.366 


1.010039766 


2 


710385528 


4 


18329885 


20 


5.797 


1.005291363 


2 


708645561 


4 


17935372 


30 


6.866 


1.003815977 


2 


708337656 


4 


17844614 


40 


7.751 


1.003106441 


2 


708230888 


4 


17805574 


50 


8.518 


1.002691899 


2 


708181518 


4 


17784097 


60 


9.202 


1.002421030 


2 


708154647 


4 


17770589 


70 


9.825 


1.002230636 


2 


708138397 


4 


17761342 


80 


10.40 


1.002089737 


2 


708127815 


4 


17754632 


90 


10.93 


1.001981392 


2 


708120532 


4 


17749550 


100 


11.43 


1.001895574 


2 


708115301 


4 


17745573 


500 


69.11 


1.000989809 


2.708093424 


4.17706229 



In Figure [T] we have plotted the wave function of the ground state of the fractional harmonic oscillator for different 
values of a and using a grid with N = 50. The case a = 2 is the exact Gaussian wave function of the standard 
harmonic oscillator. 

In Figure [2] we have plotted the energies of the first two states of the fractional harmonic oscillator for different 
values of a, using a grid with N = 60. It also shows the results given by the WKB formula derived by Laskin 0]. 
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FIG. 1: Ground state of the fractional harmonic oscillator for different values of a, obtained with N = 50. 




FIG. 2: Energies of the first two states of the fractional harmonic oscillator as functions of a. Dots and lines mark present 
results for N = 60 and the WKB ones, respectively. 



The agreement for the first excited state is remarkable, and the large deviation for the ground state is not surprising 
because the WKB method is expected to be valid for sufficiently large quantum numbers. 

Finally, we show results for the anharmonic oscillator V(x) = x A with a fractional Laplacian corresponding to a = 
4/3. We choose this particular example because the WKB formula predicts its energy spectrum to be evenly spaced, 
like the standard harmonic oscillator. Figures [3] and |4] display the energies of this fractional anharmonic oscillator 
and the absolute value of the first ten wave functions. The spectrum follows the straight line E n = 0.941 + 1.886 n in 
good agreement with the WKB estimate 

2tt / 1\ 2tt / 1\ / 1\ 

E n = ——, ; r(n+ - = : , N , rr I n + - « 1.88562 x (n+-). 

51/4,3 4+ 1) V 2 r(l/4 r 7/4 V 2 \ 2 
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FIG. 3: Energies of the fractional quartic anharmonic oscillator with a = 4/3, calculated using LSF2 with N = 50 (points). 
The solid line is the least-squares fit E n = 0.941 + 1.886 n. 





20 



































































FIG. 4: Absolute values of the wave functions (absolute value) of the fractional quartic anharmonic oscillator with a — 4/3 
and p = 4 calculated using LSF2 with N — 50. 



B. Fractional Mathieu equation 



In order to illustrate the application of the present collocation method to a problem with other boundary conditions 
than those discussed earlier, we consider the fractional extension to the well known Mathieu equation [17j | 



fv 

dz 



•j| + (o- 2qcos2z)y = 0, 



given by 



(-A) a/2 y - (a - 2qcos2z)y = 0, 



(15) 



(16) 
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which reduces to (fT5|) when a = 2. We resort to the set of periodic boundary conditions that apply to the little sine 
functions LSFi. In this case L is fixed to L = it. 

We follow the standard notation and denote by ai, i = 0, 1, 2, . . . [13], the eigenvalues with even periodic functions, 
and by i = 1,2,..., the eigenvalues corresponding to odd periodic functions. Note that, in the case of q = 0, ao = 
and the eigenfunction is constant. 

Table [TT| shows the convergence rate of our method for the lowest eigenvalue ao of the fractional Mathieu equation 
with a = 1, a = 3/2, and q = 1. The rate of convergence for the Mathieu equation is considerably larger than the 
one discussed above for the fractional harmonic oscillator. In both cases it increases as a approaches the ordinary 
value a — 2. 



TABLE II: Rate of convergence for the eigenvalue (ao) of the fractional Mathieu equation. 



N 


a = 1 


a = 3/2 


10 
20 
30 
40 


-0.78002010749909950036806303597771247593586508305415 
-0.78002010679715466707531498072556654488485913105173 
-0.78002010679715466707518738897487326807118147127619 
-0.78002010679715466707518738897487326774021967916894 


-0.60337681905510495225302969295010866313871477789813 
-0.60337681905490085108768066745917280899338056699081 
-0.60337681905490085108768066745913456676371052157098 
-0.60337681905490085108768066745913456676371052157098 


exact 


-0.78002010679715466707518738897487326774021967916894 


-0.60337681905490085108768066745913456676371052157098 



Table HTT1 shows the four lowest eigenvalues of the fractional Mathieu equation for different values of a, again with 
= 1, as well as the symmetry and periodicity of each wavefunction. They agree with the known results when a = 2 
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TABLE III: Eigenvalues ao, &i, ai, 62 of the fractional Mathieu equation with q — 1 and different values of a (N = 50). 



a 


ao 


6l 


ai 


02 


1 

3 
2 

2 

5 
2 

3 


-0.78002010679715466708 
-0.60337681905490085109 
-0.45513860410741354823 
-0.33549116582363455500 
-0.24308662756250760871 


-0.31981501215423234713 
-0.18880108186701679596 
-0.11024881699209516991 
-0.06396091681659914089 
-0.03699729990815279808 


1.2959422293970261239 
1.7046089276653617549 
1.8591080725143634723 
1.9267035413113906794 
1 .9600508496994480694 


1.5491290256879243036 
2.6389530962188063857 
3.9170247729984711867 
5.6189308675791269007 
7.9821470161415594702 


Symmetry 


even 


odd 


even 


odd 


Period 


2 7T 


7T 


2 7T 


7T 



Figure[5]shows the eigenvalues di,i = 0, 1, 2, 3, bi, i = 1,2, 3, for the fractional Mathieu equation with a = 2, 3/2, 5/2 
obtained with the periodic set LSFi and N = 30 for q £ [0, 15] (cf. [13, Figure 20.1]). Observe that the overall pattern 
of the behavior of these eigenvalues with q is similar for all values of a. This fact is hardly surprising because the 
degeneracy at q = and q — > 00 is a consequence of the form of the periodic potential. 

In this extremely simple case it may be more practical to resort to the standard method [Tt} because the calculation 
of the matrix elements of cos 2z does not offer any difficulty. We have chosen this example simply as a test; however, 
one will notice the advantage of the collocation method in the case of an arbitrary periodic potential V(z) where the 
calculation of the matrix elements may not be so simple. 



V. CONCLUSIONS 



We have devised a numerical scheme based on collocation which allows one to solve the fractional Schrodingcr 
equation on a uniform grid. We have applied this method to obtain accurate energies and wave functions of a 
fractional harmonic oscillator and compared the former with the WKB ones 0]. Our results confirm that the WKB 
approach yields reasonable results for excited states. In the case of the fractional harmonic oscillator the WKB formula 
does not give the exact result for all quantum numbers as in the ordinary case a = 2. 

We have also studied a fractional anharmonic oscillator for which the WKB formula predicts equally spaced levels 
like the standard harmonic oscillator. Our accurate results confirm this prediction beyond any doubt. Finally we have 
solved a fractional Mathieu equation, with periodic boundary conditions, and obtained the eigenvalues for different 
values of a and potential strength. 

One of the main advantages of the collocation methods in general is that they bypass the problem of calculating 
the matrix elements of the potential. This feature is most welcome when the potential-energy function is rather 
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FIG. 5: Eigenvalues of the fractional Mathieu equation for three values of a. 



complicated. This is not the case of the models chosen here because we have been mainly interested in the discussion 
of the fractional kinetic energy. 
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